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ABSTRACT 



Context. We report on the results of calibrating and simulating the instrumental polarization properties of the ESO VLT adaptive 
optics camera system NAOS/CONICA (NACO) in the Ks-band. 

Aims. Our goal is to understand the influence of systematic calibration effects on the time-resolved polarimetric observations of the 
infrared counterpart of the Galactic center super-massive black hole at the position of Sagittarius A* (Sgr A*). 
Methods. We use the Stokes/Mueller formalism for metallic reflections to describe the instrumental polarization. The model is com- 
pared to standard-star observations and time-resolved observations of bright sources in the Galactic center. The differences between 
calibration methods are simulated and tested for three polarimetric Ks-band light curves of Sgr A*. 

Results. We find the instrumental polarization to be highly dependent on the pointing position of the telescope and about 4% at 
maximum. Given the statistical uncertainties in the data acquisition, the systematic effects of the employed calibration method are 
negligible at high-time resolution, as it is necessary and achieved for in the case of Sgr A*. We report a polarization angle offset of 
13.2° due to a position angle offset of the /1/2-wave plate with respect to the header value that affects the calibration of NACO data 
taken before autumn 2009. 

Conclusions. With the new model of the instrumental polarization of NACO it is possible to measure the polarization with an accu- 
racy of 1% in polarization degree. The uncertainty of the polarization angle is < 5° for polarization degrees > 4%. For highly sampled 
polarimetric time series we find that the improved understanding of the polarization properties gives results that are fully consistent 
with the previously used method to derive the polarization. The small difference between the derived and the previously employed 
polarization calibration is well within the statistical uncertainties of the measurements, and for Sgr A* they do not affect the results 
from our relativistic modeling of the accretion process. 

Key words. Instrumentation: polarimeters, techniques: polarimetric, Polarization, infrared: general, black hole physics, Galaxy: 
center 



1. Introduction 

The polarization of electromagnetic radiation is an essential 
piece of information to determine the nature of emission pro- 
cesses and the physical parameters of the environment in which 
the radiation is generated. This is true in particular for the emis- 
sion from Sagittarius A* (Sgr A*) at the center of the Milky 
Way. This source is associated with the nearest super-massive 
black hole candidate (~ 4 x 10 6 M o ), as inferred from motions 
of stars near the Galactic center ( Eckart & Genzellll996l [l997; 
Eckart et alJl2002t ISchodel et al.H2002t Eisenhauer et all £003: 
Ghez et al.ll2000ll2005Ll2008l:lGiUessen et al.ll2009l) ~ 

Since the first near-infrared (NI R) polarimetric Wo llaston 
prism observation of Sgr A* in 2004 (Eckar tet al.l2006a[). pola r- 
ized flares have been regularly observ ed (Meve ret alj|2006allbl 
I2007t Eckart etalj|2008al IZamaninasab et al.l |201Q|) . The NIR 
counterpart to Sgr A* is extremely variable with short bursts of 



increased emission exceeding 5 mJy, which occur four to six 
times a day and last typically for about 100 minutes. The linear 
polarization degrees can reach 20% to 50% of the total intensity. 
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2008; lDodds-Eden et al.l 2009; Sa bha et alj|2010l) . This strongly 
suggests synchrotron-self-Compton (SSC) or inve rse Compton 
emission as the responsible radiation mechanis m (Eckart et al.l 
I200H l2006ald : lYuan et al.ll2004l; Eiu et alj|2006l) . 

Some models that have been applied successfully to the 
observations assume the flare phenomenon to be linked to 
emission from single or multiple hot spots near the last 
stable orbit of the black hole. The characteristic behav- 
ior of general relativistic flux modulations that are pro- 
duced v ia such orbiting hot spots have been discussed earlier 
(see e.g. lCunningham & Bardeenlll973blStark & Connorsll 19771: 
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Abramowicz et al. 1 99 ll 
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. 2008. Zamaninasab et al. 20101). 



Based on relativistic models Zama ninasab et al] (12010) find a 
correlation between the modulations of the observed flux density 
light curves and changes in polarimetric data. The authors also 
confirm that this correlation is predicted by the hot spot model. 
Correlations between intensity and polarimetric parameters of 
the observed light curves and a comparison of predicted and ob- 
served light curve features through a pattern recognition algo- 
rithm result in the detection of a signature possibly associated 
with orbiting matter under the influence of strong gravity. This 
pattern is found to be statistically significant against randomly 
polarized red noise. 

The investigation of the emission from Sgr A* and the ap- 
plication of the model calculations to the observed light curves 
of polarized light decisively depend on the quality of the po- 
larization calibration. The crucial NIR polarization data were 
obtained through NACO, a differential polarimetric imager, at 
the ESO VLT UT4. This system is mounted at a Nasmyth fo- 
cus of the altitude azimuth mounted UT4 telescope, which com- 
plicates the exact polarimetric calibration. Therefore we carried 
out a detailed analysis of the instrumental properties of this sys- 
tem, determined the systematic instrumental uncertainties and 
discuss their influence on the Sgr A* measurements and the con- 
sequences for the astrophysical interpretation of the light curves 
obtained in polarized light. 

The goal of this paper is to investigate the instrumental po- 
larization (IP) mainly on a base of scientific data, which are the 
outcome of eight years of observations of Sgr A*. These data are 
optimized for the astrophysical time series analysis of Sgr A* 
with high-time resolution. They nevertheless provide enough in- 
formation to tackle two aspects of the systematic instrumental 
effects: On one hand the description of the IP and its behavior in 
absolute values with an accuracy of about 1 % in linear polariza- 
tion degree, and on the other hand the systematic uncertainties 
of the time variability of the polarimetric parameters of Sgr A*. 

The polarimetric mode of NACO considered here does not 
provide information on circular polarization, and the flat-field 
calibration data are not optimized for polarimetric measure- 
ments. Both facts necessitate complicated procedures for cali- 
brating the existing data. We describe these procedures here in 
detail. We also give ideas on how to improve the polarimetric 
calibration for future observations that do not require a high-time 
sampling, but aim to measure polarized emission to an accuracy 
of even a few tenth of a percent in linear polarization degree. 

In section [2] we give a detailed description of the instrumen- 
tal polarization of NACO in Ks-band and a correction algorithm 
using a model based on material constants of standard coatings. 
Such a model has not yet been available to the comunity for 
this telescope/camera combination and therefore could not be 
applied to observations of Sgr A* until now. 

In section [3] we compare the model with observations of 
standard stars and time-resolved polarimetric measurements of 
bright stars in the Galactic center and analyze to what degree we 
can correct the systematics. 

In section|4]commonly used calibration procedures are com- 
pared with the new calibration method. Systematic differences 
between the methods and their influence on the variability mea- 
surements of the polarization degree and angle of Sgr A* are 
discussed. 

In section [5] we summarize our results and their impact on 
the interpretation of polarimetric time series of Sgr A*. 



2. A model for the instrumental polarization of 
NACO 

NACO is an adaptive optics imager for the NIR at the Very Large 
Telescope (VL T) that is run by ESQ. For a detailed descr iption 
of NACO see |Lenzen et all (120031) . iRousset et ail (120031) . and 
Ageorges et al. (2007). It provides a mode for polarimetric dif- 
ferential imaging combining a Wollaston prism (in the follow- 
ing referred to as Wollaston), which allows for measuring two 
orthogonal angles simultaneously, and a A/2 wave plate (HWP). 
In observation periods before 2008 a wire grid mode was avail- 
able as well. To determine the instrumental polarization (IP) of 
this instrument we use the Stokes and Mueller calculus. It allows 
us to describe the influence of optical elements on the polariza- 
tion. The model presented in this section enables us to determine 
the IP as a function of the parallactic angle. It is developed fol- 
lowing a model for the IP of the Tel escopio Nazionale Galileo 
(TNG) presented in lGiro et al.l (1200 3l) . Since it is crucial for the 
following to have a clear definition of conventions and variables 
we first introduce some basic formalism. 



2.1. Polarimetry with NACO: basics and conventions 

The projection of the electric field vector E of fully linearly po- 
larized light onto a preferential direction that makes an angle i/r 
with the E-vector is given by 



= E cos(^) 



(1) 



The energy carried by the projected electromagnetic wave and 
thus the intensity is proportional to E 2 : 



I~El=Elcos 2 m 



(2) 



For partially linearly polarized light - defined by I tol the total 
intensity, P the degree of linear polarization, and qb the polariza- 
tion angle - the dependency of the intensity on the angle of the 
preferential direction is given by 



IpA ) = ^ + cos(2 [0 - cf>]) = 1(0) . 



(3) 



is measured with respect to the polarization angle reference 
that defines <p and commonly is given by the north-south axis on 
the skjfl. 1(0) is the intensity we would measure with an analyzer 
at the angle position and a transmittance of unity. Please see 
Fig. EH 

A convenient tool to describe partial polarization of incoher- 
ent light is the Stokes formalism. The normalized Stokes vector 
for partial (linear) polarization is defined as 



S = 



with 



Q = 



u 



ho, = i ^ 

Q 
u 
v 



7(0°) - 7(90°) 
7(0°) + 7(90°) 
7(135°) -7(45°) 
7(135°) + 7(45°) 



(4) 



(5) 



1 It is the direction of the linear polarization if U = (see Eq.|5) and 
also defines the orientation of the Mueller matrices below. 



2 



G. Witzel et al.: The instrumental polarization of NACO 



0.95 - 



£ 0.9 ■ 



a 0.85 - 




15 30 45 60 75 90 
angle of incidence [deg] 



15 30 45 60 75 90 
angle of incidence [deg] 



Fig. 1. Reflection properties of metallic surfaces in Ks-band. The 
green curves represent gold, the blue aluminum. 



Q and U represent the linear polarization and V the circular po- 
larization (NACO does not provide a /l/4-wave plate to measure 
V). The parameters P and <p are related to Q and U by 



p = ve 2 + u 2 



1 (U 
— arc tan — 

2 \Q 



(6) 



In this formalism the influence of any optical element on the 
intensity and polarization can be expressed by a linear operation 
on the Stokes vector: 



S' = MxS 



(7) 



where M is the Mueller matrix of the element, S the Stokes vec- 
tor of the incoming light and S' the Stokes vector of the outgoing 
light. The elements of the Mueller matrix represent the linear de- 
pendency of each Stokes parameter in S ' on those in S . 

The optical elements considered here are mainly mirrors 
with metallic coatings. For these surfaces every incident elec- 
tromagnetic wave can be decomposed in a component parallel 
to the plane of incidence and one perpend icula r to thi s plane . 
As described in detail bv iGiro et all d2003h and lClarkd (Il973l) . 
reflections at metallic surfaces have two effects: (1) a reflection 
introduces a linear polarization because the reflectivity for these 
components is different; (2) the reflection causes a circular polar- 
ized contribution by shifting the phase between the components. 
Both effects can be described by a Mueller matrix that combines 
the matrix elements for a linear polarizer and a retarder: 



R 



' kr± + m) kr± - m) 
2</-l - n) 2^ + Hi) 










/r ± q cos((5) -^JFJ)\sm(6) 

TFlHi sin(£) yjr±T\\ cos(5) 



(8) 



with r ± and r\\ the reflection coefficients for the two components 
and 6 the relative retardation between the components. With the 
material-dependent refractive index and extinction coefficient all 
three parameters can be calculated by using the Fresnel formu- 
lae. In the left plot of Fig. [T] we show the reflection coefficients 
for p- and s-waves as a function of the angle of incidence for Ks- 
band. In the right plot of Fig. [T] we show the relative retardation 
a mirror hit by a perpendicular beam is considered to have the 




detector 



II ® 




)J2 plate 



SI 3/27 optics 



Fig. 2. Optical elements of UT4, NAOS & CONICA and their 
relative orientation in the moment of the meridian transit. 



same effect as a HWP. The material constants (Ks-band) used for 
these plots are listed in Table [2] 

Eq. |8] defines the Mueller matrix for metallic reflection in a 
way that the preferential direction of the matrix (the direction of 
the introduced polarization, always perpendicular to the plane of 
incidence) is oriented parallel to the polarization angle reference 
(generally north-south). To change the (perpendicular) orienta- 
tion of the plane of incidence with respect to the polarization 
angle reference by an angle y, one has to apply the transforma- 
tion 



R' = T(-y) xRx T(y) , 

with T the rotation in Stokes space, which is defined as 



T(P) 



1 0^ 

cos(2p) sin(2/?) 

~sin(2/>) cos(2/?) 
1 ) 



(9) 



(10) 



NACO is mounted at the Nasmyth focus of VLT Yepun. 
Because of the 45° tilted folding mirror M3 (see Fig. |2| that 
sends out the beam to the Nasmyth focus, NACO has a signifi- 
cant instrumental polarization that is depending on the parallac- 
tic angle. In the ESO user manual the total IP is estimated to be 
up to 4%. NACO provides various setups for polarimetry com- 
bining different cameras (with different pixel scales) and filters 
with the Wollaston an d the HWP. A table o f available cameras 
and filters is shown in Ageo rges et alJ (120071) . Here we are con- 
centrating on the most commonly used setup: The cameras S13 
and S27 with the Ks-band filter, the Wollaston, and the HWP. 

The field of view of S 13 in combination with the Wollaston 
analyzer and the polarimetric mask is shown in Fig. [3] The upper 
stripe is the ordinary beam (0°) and the the lower one the extraor- 
dinary (90°). The angle positions of the HWP in the headeiQ of 
the NACO data is counted with the same sense of rotation. Note 
that this sense is reverse to the sky because the number of mirrors 
is odd and every mirror turns it once (see Fig.|2]and Fig.[4]l. 

Not only M3 contributes to the IP of NACO. Every signifi- 
cantly inclined reflective surface in the light train has to be in- 
cluded in a model of polarimetric instrumental systematics. All 
optical elements (including the analyzer) and their relative ori- 
entations are discussed below. In Fig. [2] we show the optical el- 
ements of UT4 and NACO. In this sketch the direction of the 
optical train within NAOS (defined by the connection line be- 
tween the input mirror and PI in Fig. |4|i is perpendicular to the 



2 The angle position of the HWP is reported in the NACO FITS 
header under the keyword "INS RETA2 ROT"; the encoder position 
can be found under "INS ADC1 ENC". 
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PI 



Fig. 3. Polarimetric observations with the Wollaston: the picture 
shows a frame observed with the S13 camera. The image taken 
with the ordinary beam of the Wollaston is the upper stripe. 



paper plane and indicated by the face-on arrow sign on NAOS. 
A sketch of NAOS is shown in more detail in Fig. [4] 

A de-rotation of the parallactic angle rotates the optical train 
beginning with NAOS with respect to the mirrors Ml/2/3. The 
elevation rotation is only affecting the hole assembly because the 
rotator of NACO compensates for it. It does not influence the IP. 
After HWP and Wollaston there are two more folding mirrors in 
the architecture of CONICA; the inclinations of these two mir- 
rors are different for different cameras (see Table [2] and angles 
of incidence e therein). The light train of NAOS in Fig.|4]shows 
that there are only two mirrors with a significant inclination: the 
input and output mirror (red). The parabolic mirrors PI and P2, 
the tip-and-tilt mirror TTM, and the deformable DM have in- 
clinations < 5° and can be neglected. The specifications of the 
dichroic are not accessible to us, but here the inclination of only 
about 12° is also very small. 

2.2. Instrumental polarization generated by M3 

M3 is coated with aluminum and inclined by 45°. To understand 
the position-dependent part of the IP, it is important to analyze 
the time-depending orientation of M3. 

The common angle reference at the sky is the north-south 
axis. At the Nasmyth focus this direction is parallel to the plane 
defined by the edge of the main mirror Ml in the moment of the 
meridian transit of the source. In this moment the north-south 
axis is parallel to the preferential direction of the matrix in Eq.[8] 
(see Fig. 0. At all other moments it is tilted by —p with respect 
to the preferential direction of matrix (O (with p the parallactic 
angle). 

In order to transform matrix (O to the north-south reference 
at any moment, we have to use a transform similar to Eq.|9] After 
changing the reference from the celestial system to the reference 
system of M3 and applying R for the metallic reflection, we have 
to transform to the detector reference. The detector is de-rotated, 
it follows the elevation rotation of the Nasmyth focus and the 



NAOS 



DM 



output 



input 



dichroic 



TTM 



P2 

Fig. 4. Light train in the adaptive optics module NAOS 

parallactic rotation of the sourc^- After the reflection this paral- 
lactic rotation has the reverse sense, thus the position depending 
part of the IP can be expressed as 



M 



T{p) x R alu x T{p) 



(ID 



where R a i u is the reflection matrix for bare aluminum with the 
values listed in Table [2] 



2.3. A description of the entire instrumental polarization of 
NACO 

NAOS 

The two 45° inclined Silflex coated folding mirrors in the adap- 
tive optics module NAOS are described by the square of ma- 
trix ®: 



NAOS 



— Rril X R. 



(12) 



The material constants for Silflex have been provided by the pro- 
ducer Balzers Optics and can be found in Table [2] 



The HWP 

The polarimetric analyzer is part of the camera CONICA. It 
mainly consists of the HWP, the Wollaston, and the detector. The 
HWP turns the angle of the linearly polarized part of the light by 
the double of its position angle. The position angle is the angle 
between the fast axis (or slow axis, the degeneracy is 90°) and 
the polarimetric angle reference. The formula (see TableQ]) given 
in the ESO manual for the dependency of the position angle (an- 
gle with respect to the north-south axis) on the encoder steps has 
to be modified. The plot in Fig.[5]shows the un-calibrated polar- 
ization of IRS21 in a dataset taken in 2005 with the Wollaston. 
The dataset exhibits an offset in pola rization angle of about 14° 
in comparison with lOtt et al.l (1 19991) . A maintenance of NACO 



i.e. the angle between main mirror and detector orientation only 
depends on the parallactic rotation. 
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Table 2. Material constants of coatings. 
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Fig. 5. IRS21 in 2005 (Wollaston). The plot shows the intensity 
divided by I, or /2 as a function of angle 8 (as defined in Eq. [3]). 
The data are obtained from mosaics. No calibration was applied. 
The green line marks the measured po larization angle, the blue 
the published angle (9.8%@14°) from lOtt et al] dl999l) . In this 
plot the angle is counted in the instrument sense (negative with 
respect to sky). 



Coating Mirror n 


k 


e 


11 


' perp 


6 






[deg] 






[deg] 


Alu. M3 2.75 


20.0 


45° 


0.96262 


0.98113 


176.03 


Silflex NAOS I&II - 




45.00 


0.98272 


0.98872 


165.00 


Gold S13 I 0.99 


13.8 


42.05 


0.97266 


0.98484 


175.02 


S13 II 




2.95 


0.97959 


0.97964 


179.98 


S27 I 




27.65 


0.97701 


0.98193 


177.99 


S27 II 




17.35 


0.97865 


0.98054 


179.23 



Notes. Refractive index n, extinction coefficient k, reflection coeffi- 
cients, and relative retardation for the coatings of NACO at 2.2//. e is 
the angle of incidence for which the reflection coefficients and rela- 
tive retardation was computed. The material constants can be found on 
www.RefractiveIndex.info. The default k value for aluminum is 22.3, 
and <5 for Silflex is 166.6 according to Balzers Optics specifications. 
Both values have been changed to match the position dependency of 
measured Stokes parameters (see section [3~2l . 



Table 1. Encoder positions of the HWP. 



Encoder steps Angle 
[deg] 



Manual : 






X 


a = (x + 205) * 


0.08789 


-205 ~ 3891 







51 


22.5° 







205 * 0.08789 = 


= 18.02 


Revision : 









(11.4 ±0.2) 




Difference : 








(-6.6 ±0.2) 





Notes. The table shows that the actual reference system for the HWP 
is offset by (-6.6 ± 0.2)° with respect to the reference system assumed 
in the manual. This results in a positive angle offset of 13.2° for the 
polarization channels (in the instrumental sense of rotation). One turn 
of the plate corresponds to 4096 encoder steps. An angle of 0° actually 
corresponds to an encoder position of 3966. 



Total IP 

Because the plane of incidence of the NAOS light train is per- 
pendicular to the plane of incidence of M3, we have to transform 
Mnaos according to Eq.|9]with y = 90°. We also have to add a 
rotation matrix for the rotator adapter that de-rotates the instru- 
ment and can change the orientation of the field of view (from 
north-south on the y-axis of the detector to any angle). The dis- 
cussed matrices and orientations finally result in a total Mueller 
matrix: 

M NAC0 = (Tr x) T(-/3) x 7\90°) x M NA0S 
xT(-90°) x T(a) x M M3 
= (Tr x) T(-p) x 7\90°) x R sU x R sU x 7\-90°) 
xT(a)xT(p)xR alu xT(p) , (13) 

with -a the angle of rotator adapter as reported in the ESO 
FITS header keyword "ADA POSANG" and /3 = 13.2° the 
offset of the HWP. The matrix Tr represents the effects of the 
analyzer and its transmission. It is included here for the sake of 
completeness. A more detailed discussion of these effects and 
their correction will follow in section 1231 



in autumn 2009 revealed the actual position angle reference. The 
true offset is (13.2 + 0.3)° (see TableQ]), a value that agrees very 
well with the value of 14° that we predicted from the observa- 
tional data before t he intervention. T he previously reported an- 
gular offset of 34° |Trippe et al.l2007l) could not be verified. This 
value may result from th e offset we found combined with the fact 
thatUrippe et al. (2007) use a sine- rather than a cosine-function 
to describe the 20-dependency of the polarized channel flux. This 
causes a 45°-shift in the opposite direction to the HWP offset, re- 
sulting in a total of 32°. 

The HWP was installed on 2003 August 8, and we assume 
that the offset is constant for all epochs since then. 



2.4. The instrumental polarization in numbers 

In this section we investigate the behavior of the introduced 
model. All material-dependent parameters of this model are 
summarized in Table [2] These parameters are mainly literature 
values for the materials. In section[3]the model is compared with 
standard stars and light curves of bright GC stars, which exhibit 
the variations of the IP with the parallactic angle. To match the 
observations the material constants k for aluminum and 6 for the 
Silflex coating had to be slightly changed as described in the cap- 
tions of Table [2] Here we already discuss the final version of the 
model that is gauged with the calibration sources. In section [3] 
we will then justify the model and the chosen parameters. 

We are now able to evaluate the contributions of the differ- 
ent optical elements quantitatively. First we numerically express 
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the matrices of Eq. [13] For the given material parameters and a 
parallactic angle of 0° we find Mmt, to be 



0.972 0.009 
0.009 0.972 
-0.969 -0.067 
0.067 -0.969 



(14) 



and for a parallactic angle of 45° 



M M3 = 



0.972 0.009 
0.969 

-0.009 -0.972 
-0.067 





-0.067 


-0.969 



(15) 



Obviously the main effect of the tertiary mirror is an / — > 
cross talk of about 1 % of the total intensity. The cross talks be- 
tween linear and circular polarization can be on the order of 7% 
of Q/U or V respectively. For weakly polarized sources this is a 
minor effect. The transformed M^aos is given by 



7X90°) x M NA0S x T(-90°) 



I 0.972 -0.006 

-0.006 0.972 

0.841 -0.486 

0.486 0.841 



.(16) 



Here Q is affected by about 0.6% of / and the U <-> V cross 
talks are on the order of 50% of the corresponding values U or 
V respectively. With these numerical matrices, Eq. Qj] and the 
parameter a set to zero, we can write matrix M NAC o as 



0.944 

-0.005 + 0.008c ; , 
-0.003 + 0.004c,, • 
-0.004.?,, 

0.009 s p -0.001 
0.366 + 0.015c p 
-0.731 -0.029c,, 
-0.472 + 0.057c„ 



0.009c p - 0.006 
+ 0.003j p 0.843 -0.015^ 
- 0.007 s„ 0.422 + 0.029^,, 



-0.057i-„ 



-0.211 



0.000 
• 0.059.S 



„ + 0.025c p 



0.421 -0.029.S,, -0.051c,, 
-0.816 - 0.033c„ 



(17) 



with c p - cos(p) and s p = sin(p) describing the dependency 
on the parallactic angle. Tr was not included. This last expres- 
sion gives all necessary information on the cross talks and their 
dependence on the parallactic angle. The cross talks from I to 
the linear polarization is on the order of 0.5% - 1% of the total 
intensity. There are also strong interactions between Q and U 
(on the order of 40% of the corresponding value Q or U respec- 
tively), between U and V (also on the order of 40%) and from 
V to Q (on the order of 20%). For not very strongly polarized 
sources the position-dependent variability of the linear polariza- 
tion is dominated by the I <-> Q/U cross talks. 

As an example for the telescope position-dependent behav- 
ior of the IP we plot in Fig. [6] Q, U, V, the linear and total po- 
larization degree, and the polarization angle as functions of hour 
angle for an unpolarized source at the position of the GC (May, 
alt = 85.6° for meridian transit; solid lines in the plots). The IP 
of NACO reaches about 1.6% at maximum. Around hour angle 
zero the IP changes most rapidly, as is expected. Here the polar- 
ization degree reaches its minimum and the polarization angle 
swings to its other extreme. The curves are asymmetrical around 
zero because of the HWP offset. 



4 Depending on the parallactic angle, this cross talk is affecting Q 
(p = 0°/90°), U (p = 457135°), or both of them. 




Fig. 6. Upper panel: Instrumental polarization predicted by the 
model for an unpolarized source corresponding to matrix (TTTb 
(solid line) and including the systematic effects of the the ana- 
lyzer (see section l2~5l and Fig.|7]therein, dashed line). The curves 
show the Stokes parameters Q (blue), U (red), and V (green) as a 
function of hour angle. Here we considered the most common in- 
strumental setup with Wollaston, HWP, and S13 optics. Middle: 
Linear (green) and total (red) polarization degree corresponding 
to the upper Stokes values. Lower panel: Polarization angle. 



2.5. The polarimetric analyzer 

Up to this point we discussed the influence of the optical ele- 
ments that are located in front of the analyzer. In this section we 
consider the analyzer itself and its systematic effects on polari- 
metric measurements. 

The main element of the analyzer is a polarizer. After the 
polarizer the only quantity of interest is intensity, because (for a 
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polarizer with an efficiency sufficiently close to 100%) the polar- 
ization is known to be 100% in the direction of the polarizer (see 
Eq. [5). Thus, the polarizer and the optical elements within the 
analyzer after the polarizer have to be investigated with respect 
to their relative attenuation of the different channels. 



CONICA 

After the Wollaston the optics of CONICA include two more 
folding mirrors. These mirrors have gold coatings and show dif- 
ferent inclinations for S13 and S27. Because after the Wollaston 
only the intensity of the two channels is important, it is not use- 
ful to work with matrix ©, which would describe the change in 
polarization inside the analyzer. We are interested instead in how 
the measurement of the Stokes vector of the light in front of the 
analyzer is affected by the different attenuation of the intensities 
of both channels inside the analyzer. 

For CONICA the relative attenuation factors for the orthogo- 
nal channels of the Wollaston are the same for the measurement 
of Q and U, because the arrangement of the optical elements af- 
ter the Wollaston does not depend on the position of the HWP. It 
is possible to find a Mueller matrix for this influence. 

The matrix to correct for a transmittance difference of the 
two orthogonal polarimetric channels after the Wollaston can be 
deduced in the following way: Let be I\ and I 2 the intensity of the 
ordinary and the extraordinary beam, respectively. The folding 
mirrors after the Wollaston attenuate these intensities: 



Wollaston flat-field, deg vs. 90 deg channel, 2005-07-30 



V 2 = nh 



(18) 



with r\ < 1 and r 2 < 1 the attenuation factors of the ordinary 
and the extraordinary beam, respectively. Thus, the measured not 
normalized Q (other stokes parameters U and V analog) will be 



Q' = A-1'2 
r = i[ + i 2 



(19) 



The correction matrix applied to the measured Stokes vector S ' 
has to give the vector S , which results from the intensities I\ 
and l 2 . Because each of the Stokes parameters Q, U, and V is 
determined from a pair of orthogonal Wollaston channels, all of 
them are affected in the same way, and we can expect e.g. Q to 
depend on /' and Q' only (analog for U and V): 



Q = xT+yQ' . 
This leads to the equation 
1 



y 



•4- 

V2 



(20) 



(21) 



which has the solution 

n - n 



y 



n + n 

2nr 2 



(22) 



In a similar way we find the matrix elements for /. The resulting 
matrix is the inverse matrix of the transmission matrix Tr of 
Eq.[T3land can be written as 



Tr 1 = 



2nr 2 



T + T 

T T + 

T T + 

T T +J 



(23) 




flat— field value 

standard flat— field, areas of deg vs. 90 deg channel, 2005-07—30 




1 1.1 
flat-field value 



Fig. 7. Histogram of flat-field pixel values in the regions of the 
orthogonal channels for the S13 polarimetric mask. The upper 
plot shows the histogram for a detector illuminated by the cal- 
ibration lamp through the polarizer, the lower one for a simple 
imaging twilight flat. 



with T ± = r 2 ±r\, and by inversion we find 



Tr = 2r\r 2 



T + _ 

+ tS 



r_ 



T k T - 

T 2_ T 2 
jl 



^ 




1 



■ T 2_j2 r+ ( r 2_ r 2) r + 
T 2_j-2 T+ ( T 2_ T 2_y " F+ J 



(24) 



For the Wollaston (polarization of the ordinary beam orthog- 
onal to the plane of incidence at the first CONICA mirror) the at- 
tenuation factors r\ and r 2 can be computecQ from the reflection 
coefficients for gold for the corresponding angles of incidence 
(for values see Table |2]l 



■old,! gold.II 



r 2 



-kit 

_ ygoldj eold.II 
— r ± r ± 



(25) 



Unfortunately the information on the attenuation caused by 
the two mirrors is not sufficient for a correction of the IP of 



5 Under the assumption of 100% efficiency of the Wollaston as a po- 
larizer and with the polarization of the ordinary beam perpendicular to 
the plane of incidence (first gold mirror of CONICA). 
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Fig. 8. Normalized transmission for polarized light para llel fen) 
and o rthogonal (k\ n ) to the wires of the grid (Fig. from Hodapp 
1984). Grid period 0.25 yum, equally spaced on a CaF2 substrate. 



CONICA. The characteristics of the Wollaston, in particular a 
possibly different transmission of the two beams, are not taken 
into account. These specifications are not accessible to us, but 
information on the Wollaston and its transmission is carried by 
the flat-field that is taken routinely with mask and polarizer. To 
correct for these transmission differences the relative weighting 
of the orthogonal channels in the flat-field has to be conserved 
while normalizing the flat-field. This relative weighting also con- 
tains information on the attenuation of the CONICA mirrors af- 
ter the Wollaston and additionally on the polarization of the cal- 
ibration lamp that should be zero, but in reality does contribute. 
The calibration lamp consists of a halogen bulb within a slot 
of an Ulbricht sphere, which effectively depolarizes the light 
(~ 10~ 3 % remaining). This light is then coupled in by a 45°- 
tilted gold mirror. This results in a maximum linear polarization 
of ~ 1% of the calibration light. 

A comparison between a typical imaging and a Wollaston 
flat-field is shown in Fig. [7] Obviously the flat-field with the po- 
larizer shows a significantly different distribution for the detec- 
tor areas of the two channels, while a flat-field without polarizer 
exhibits a comparatively homogeneous response in these areas. 
The ratio of the mean values of the areas that correspond to the 
extraordinary and the ordinary beam respectively is 1 .027 for the 
flat-field of 2005 (shown in Fig. |7), whereas the ratio of r 2 and 
r\ obtained with Eq.|25]is only 1.006. Thus, we recommend the 
correction by the Wollaston flat-field because the influence of 
the transmission of the Wollaston (which is also different within 
the field of view of the individual channels) to the polarization 
can be on the order of 2% in polarization degree. The dashed 
lines in Fig. [6] show the same model as given by the solid lines, 
but now including the effects of the analyzer (simulated using 
Eq.[T3l including Tr with rilr\ = 1.027). The general behavior 
of the curves is the same, but the maximum of the IP is now as 
high as 4% and the angle flips over by full 180°. 

We point out here that the polarization of the flat-field cali- 
bration light is on the same order as the remaining systematic er- 
ror of the model described here and probably is the major contri- 
bution to the deviation from the standard stars (see section [3721 . 

By using matrix ( |24| ) we assume that the attenuation factors 
are the same for the measurement of Q as for the measurement of 
U, which is the case for the Wollaston, but not for the wire grid 
mode. The transmission correction for wire grids is even more 
complex. For NACO four wire grids were available for observa- 
tions until and including 2007. The position angles of these grids 
can be found under the header keyword "INS OPTI4 NAME", 
which indicates the direction of the wires, thus for the polar- 



ization direction that is transmitted, 90° have to be added. For 
this mode one has to consider the more general case of a polar- 
izer with an efficiency r\ < 1. Fig. [8] shows the efficiency as a 
function of the wavelength. Let r\\ be the efficiency of the polar- 
ization direction perpendicular to the wires and r\ 2 the efficiency 
parallel to the wires. Then the expressions for x and y in (1221 are 
still valid if we replace r\ and r 2 with r\ and r' 2 , where 

r \ = mn- m r 2 

r' 2 = m r i- mn ■ (26) 

Here the attenuation factors are not the same for the measure- 
ments of the Stokes parameters Q and U : r\ and r 2 have to be 
calculated separately for the two position angle pairs (0°/90° 
and 45°/135°) of the wire grids according to their the orienta- 
tion with respect to the mirrors within CONICA. Because in this 
mode the flat-fields have also to be taken separately, they cannot 
be used easily to infer information on the relative transmission. 
The "boot strapping" method described in section 14.21 circum- 
vents these problems. 



2.6. A correction algorithm 

To obtain the true polarization of a source we can apply the in- 
verse of the Mueller matrix to the measured Stokes vector: 

S = Cnaco xS' = M n \ co x S' . (27) 

This is possible because every matrix in Eq.[13]is invertable and 
therefore Mnaco as well. An analytic solution for Cnaco can 
be found in appendix lAl With NACO we cannot gain informa- 
tion on the circular polarization, and we have to assume V = 0. 
However, because in the NIR the circular polarization for most 
sources can be neglected (V = 0) and the instrumental circular 
polarization of NACO can be on the order of 1 %, the assumption 
of V =0 introduces an error of this range to the other Stokes 
parameters when applying the correction matrix Cnaco- To over- 
come this problem, an iterative algorithm has to be implemented 
with the following steps: 

1 . Compute the normalized parameters Q' and U' from NACO 
data; complete the Stokes vector S' with /' = 1 (because 
the Stokes parameters are normalized by intensity) and set 
V = 0. 

2. Apply Cnaco ( as shown in appendix lAl and obtain a first 
guess on the corrected Stokes vector: 5,=o = Cnaco xS'. 

3. To initiate the z'-th step define S = (I, Q, U, V) by setting 
/ = /,-!, Q = Qi-u U = Ui- U V = 0. 



NACO 



4. Compute the numerical inverse of Cnaco (C na < 
and apply it to S : S ' = C N \ C0 x S . 

5. Replace V — in S' (defined in step 1) with V' + (com- 
puted in step 4) and get 5,- = Cnaco x 5". 

6. For the next iteration step start from 3. 

This algorithm quickly converges to a stable set of Stokes param- 
eters. After 10 iterations the differences in the obtained solutions 
approach the computational uncertainties. It replaces V — by 
a value that guarantees V to be zero. 

The algorithm has to be applied for every single frame and 
its corresponding parallactic angle and rotator position. Often it 
is necessary to work with mosaics that are obtained by averaging 
over a number of frames. Since the differences of the elements 
of the correction matrices for individual frames are often small 
within a dataset, it can be a suitable approximation to use the 
average of the matrices to correct polarimetric data that are ob- 
tained from mosaics. 
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3. Observations 

In this section we justify the presented model. The comparison 
with calibration star data shows that we can indeed describe the 
parallactic angle-dependent IP as it is observed. It also reveals to 
which accuracy we are able to compensate the systematic effects 
with our model. 

In principle it would be possible to determine the instrumen- 
tal polarization of an optical train by calibration measurements 
only: sufficiently bright sources with zero polarization, with lin- 
ear polarization at 0° (U = 0), with linear polarization of an 
angle close to 45° (U + 0), and with significant circular polariza- 
tion (or simply four linearly independent Stokes space elements) 
would provide the necessary information to solve the equations 
for the matrix elements. In practice it is not feasible to mea- 
sure four sources for every parallactic angle, rotator position, 
and both optics, and in particular it is not possible to measure 
the Stokes parameter V. However, the results of the presented 
model can be compared with the small number of available and 
suited standard observations in the ESO archive. In this way it is 
much easier to obtain a calibration of the IP. 

3.1. The data 

In order to gauge the model for the instrumental polarization we 
have to test the predictions of the model against standard star ob- 
servations. In the data archive of ESO plenty of standard star ob- 
servations are available, mainly of unpolarized standards. Most 
of these data are difficult to handle because of bad weather condi- 
tions, insufficient brightness, or the small number of frames that 
have been observed. Furthermore polarimetric measurements at 
different angles have been obtained for these standards by turn- 
ing the whole instrument instead of the HWP, which changed 
the IP. Thus, in these cases Q and U have been measured with 
different instrumental setups. However, three unpolarized stars 
turned out to be suited to be compared with the prediction of our 
model. A fourth, the standard RCra88, was even observed with 
the sa me observing s trateg y as Sgr A*. This star has been stud- 
ied bv lWhittet etaD dl992h . and it shows a K-band polarization 
of (1.8 ±0.1)% at (95 ±1)°. 

The available polarimetric standard star data in the archive 
do not include longer light curves, which would allow us to in- 
vestigate the position dependency of the instrumental polariza- 
tion. For this purpose we use bright stars in the IRS 16 cluster of 
the Galactic center that were observed with a time sampling of 
about 4 min for Stokes parameters I, Q, and U. 

In order to finally test the systematic effects of different cal- 
ibration methods on the polarimetric light curves of Sgr A*, we 
investigated three of the brightest polarized flares in the frame- 
work of the new method. All data sets used for this paper were 
taken in the Ks-band and are listed in Table IB. ll 

For the GC observations the infrared wavefront sensor of 
NAOS was used to lock the AO loop on the NIR bright (K-band 
magnitude ~ 6.5) supergiant IRS 7, located about 5.6" north of 
Sgr A*. For the standard RCra88 the AO was locked on the tar- 
get itself. For all other standards the AO loop was open. 

During the observations the atmospheric conditions (and 
consequently the AO correction when the loop was closed) were 
stable enough for high angular resolution photometry and po- 
larimetry (typical coherence time > 2ms). The exposures were 
jittered by a few arcseconds. All frames were sky-subtracted, 
flat-fielded, corrected for bad pixels, and aligned wit h sub -pixel 
accuracy by a cross-correlation method dDevilrardlll999l) . For 
the GC observations point spread functions (PSFs) were ex- 



Table 3. Standard stars. 



Source 


P 




|model - data| 




[%] 


[deg] 


[%] 


RCra88 


1.8 


95 


(0.2 + 0.2) (Q&U) 


WD 1344 







(0.3 ± 0.2) (g) 


WD2039 - 202 







(0.3 ± 0.2) (Q) 


HD 109055 







(0.3 ± 0.2) (0 



Notes. Reference polarizations of the standard stars and systematic de- 
viation of the data from the model. For RCra88 the data allowed us to 
compare both parameters Q and U. For all other standards only the dif- 
ferential flux of the orthogonal channels (equivalent to Q) for different 
instrumental setups could be tested. 

tracte d from the individual frames with StarFinder (Diolaiti et al. 
2000), the images were de-convolved with the Lucy-Richardson 
(LR) algorithm (which is necessary to counter source confusion 
and crowding), and a beam restoration was carried out with a 
Gaussian beam of a FWHM corresponding to the resolution at 
2.2yura (~ 60 mas). 

As a preparation for the differential measurements, all chan- 
nels were aligned to each other on subpixel-scale. Flux densities 
were measured by aperture photometry. The radius of the aper- 
tures for the de-convolved GC data was about 40 mas (3 pixels, 
S 1 3), about 270 mas (20 pixels, S 1 3) for the isolated RCra88 and 
about 600 mas (22 pixels, S27) for all other standards. Because 
of the channel alignment the positions of the apertures were 
the same for all channels. As a correction for background flux 
we subtracted the flux measured in apertures where no appar- 
ent source is located for each channel. Total intensity for the 
light curves was obtained by adding the flux of the orthogonal 
channel and applying a fl ux density calibration as described in 
IZamaninasab et all (l2010» . The polarimetric parameters Q and 
U were then obtained according to Eq. [5] 

3.2. Gauging the model: standards and IRS16 stars 

The GC observations and the measurement of RCra88 used the 
Wollaston-HWP setup, whereas for all standard-star observa- 
tions NACO itself was rotated for the different angles. In the 
former case it was possible to compare the data (after correcting 
for the IP of CONICA by the flat-field as described inl2~3land 
for the offset of the HWP) in Q and U with the predictions of 
the model. In the latter case the IP was different for the different 
angles of the instrument rotation. Thus, a measurement of Q and 
U under the same circumstances is not available. In these cases 
we computed the expected difference between the fluxes of both 
channels for each frame and its instrumental setup and compared 
it with the data (i.e. we considered each frame as a measurement 
of Q with a different instrumental setup). For all polarimetric 
standards the model agrees with the observations with an accu- 
racy of below 0.5 % (in Q and Q&U respectively). The standard 
stars, their reference polarizations (degree and angle), and the 
systematic deviation of the model from the differential fluxes in 
the observations are given in Table[3] The uncertainty of this sys- 
tematic error was determined taking the median deviation of the 
data, i.e. the statistical error of the measurements. 

The description of the position-dependent part of the IP can 
be tested with the long light curves of bright IRS 16 sources from 
2009. Here a comparison revealed that k" lu and 8 had to be 
slightly adjusted to match the shape of the observed light curves 
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in Q and U as described in the captions of Fig. [2] This excursion 
from the default material constants is within the typical toler- 
ances and can probably be explained with an aging of the alu- 
minum coating. The time-dependency of the model and the data 
for Stokes Q and U excellently agree with each other. Since the 
apparent polarization of sources at the GC is not known with 
the accuracy of standards, the apparent polarization parameters 
of these sources are considered to be the free parameters of the 
fits. 

Fig. [9] shows model and data for IRS16C (In Figures IbTTI 
we show the same model vs. data comparison for IRS16NW, 
IRS16CC, and S67). In this figure Q (blue) and U (green) pa- 
rameters are shown as a function of hour angle. The solid line 
is the best ^ 2 -fit of the model to the 2009 data that is already 
corrected for the HWP offset (upper panel). The lower panel 
shows the data corrected for the full instrumental polarization. 
The solid line here describes Q and U corresponding to the ap- 
parent polarization of the source. The big difference between the 
measured and the corrected U parameter results mainly from the 
opposite sense of rotation of the HWP with respect to the sky. 
The larger deviation of the model from the data points at the be- 
ginning of the night are due to the weather conditions and the 
AO performance, which were not as stable as later. The fitting of 
the model is weighted toward the end of the night. 

For all four fitted sources the errors of the single data points 
are about 0.3% for Q and U and^ 2 /dof = 1.1, where p and 
(p (apparent polarization), and k" lu and 6 s ' 1 (material constants, 
see Table |2]l are considered to be the free parameters. The ap- 
parent polarizations as best^ 2 fitting results are listed in Table|4] 
for all four IRS 16 sources obtained from 2009 data. In this ta- 
ble the error of the polarization degree is about 0.8%, the er- 
ror of the angle about 3°. For comparison we list the results of 
lOtt et alJ dl999t) . Here p has in average an error of 2%, and the 
angle uncertainty is about 18°. Additionally we give the aver- 
age polarization of the central arcsecond around Sgr A* for stars 
with rrik > 13. Stars of this brightness in the near surrounding of 
Sgr A* have commonly be used as calibration stars for the "boot 
strapping" calibration described in section 14.21 We emphasize 
that to our knowledge the results presented here are the first po- 
larimetric measurements of source s within the central pc of the 
GC since Kna cke & Cappsl (1 19771) that are independently cali- 
brated with a method that goes beyond "boot strapping" proce- 
dures (see section l4~2l . The polarization of the sources in Ta ble[H 
compa res well with the polarization found by Knacke & Capps 
dl977l) for the central region of the GC and can therefore be ex- 
plained by the galactic fore ground polariza tion as discussed in 
iKnacke & Cappsl (1 19771) and lott et alJ (Il999h . 

The presented model enables us to correct the IP with an ac- 
curacy better than 1% in polarization degree and better than 5° 
in polarization angle for polarization degrees > 4%. These er- 
rors are deduced from the light curves of the IRS 16 sources that 
exhibit comparably small statistical errors for Q and U, and the 
systematical deviation of Q and U of 0.4% for the standard stars 
(see Table [3]). The time variability (i.e. the relative behavior) of 
the IP in polarization degree can be described with an accuracy 
of a few tenths of a percent. Thus, after correction of the IP, re- 
maining variability with amplitudes of 1 % or more in linear po- 
larization is caused by intrinsic variability and statistical errors 
and is not a feature of the IP. 



6 By "apparent" polarization here and below we mean the polariza- 
tion that an ideal instrument would measure. The polarization of the 
stars discussed here is dominated by the interstellar medium. 



Table 4. Apparent polarization of sources at the GC. 



Source 


This 


paper 


Ottetal. (1999) 






P 


4> 


P 


<P 






[%] 


[cleg] 


[%] 


[deg] 




IRS 16C (S96) 


4.6 


17.8 


4.0 


35 


9.55 


IRS16CC 


5.4 


15.5 


6.1 


54 


10.15 


IRS 16NW (S95) 


5.9 


12.0 


4.6 


24 


9.86 


S67 


5.2 


17.8 






12.10 


central arc sec 


5.3 


27 









Note s. K-band mag nitudes are taken from lOtt et al.l ( 1999) and 
lGillessenetai1ll2009]) . 




-2-101234 



hour angle [ho] 

Fig. 9. Q (blue) and U (green) parameters as a function of time 
for IRS16C. The lower panel shows the data corrected for the 
full instrumental polarization. 



4. Comparison of different calibration methods 

In this section we describe two other common calibration meth- 
ods. We investigate their ability to correct the IP of NACO. Both 
calibration methods have shortcomings, and we give a quantita- 
tive analysis of the remaining systematic errors. We investigate 
the influence of the "boot strapping" method that was commonly 
used for calibrating light curves of Sgr A* on the evaluation of 
time series of the polarimetric parameters. But first we discuss 
the "channel switch" method, which enables us to propose an 
observing strategy that aims at high-accuracy polarimetry, but 
not at high-time resolution. 
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4.1. The "channel switch" method 

A common method to treat instrumental po- 
larization is to observe orthogonal channels at 
0°/90°, 457135°, 9070°, 135 745°. We assume that the IP 
affects each of the orthogonal channels factorially, and that the 
particular factor does not change if the instrument is turned by 
90°: 



£o° = M(0°) - 0-/(90°) and Q 90 ° = ,11(90°) - crl(0°) 



(28) 



with A and cr the factors introduced by the IP (for U analog equa- 
tions). The true Stokes parameters then can be obtained from the 
difference of the Stokes parameters derived from the correspond- 
ing images 0°/90° and 90° /0°, and 45°/135° and 135°/45° re- 
spectively, and with Eq.[5]we get 



Qo 



<2o° - Q 



'ro/,0 



ioi,90 



(A + cr)(/(0°)-/(90°)) 
(A + cr)(/(0°) + 7(90°)) 



Q 



(29) 



(for U analog). 

The assumption of equal factors in Eq.[28]is not appropriate 
for every case. Obviously the influence of a transmission matrix 
like Tr of Eq.|24]can be fully corrected by "channel switching", 
whereas for the matrices (0 and ( fTOb systematical errors remain. 
To quantify the remaining effects for a given Mueller matrix M 
we calculate the Stokes vectors 



So- = M x S and S' 90 , =Mx T(90°) x S 



(30) 



with S the Stokes vector of the incoming light and SL and SL. 
the measured Stokes vectors of the image pairs 0°/90°,45°/135° 
and 90° /0°, 135°/45° respectively. For matrix © we find 



Qcorr = Q and U COIT = 



2(r ± +r,|) 



-U , 



(31) 



with Q and U the proper Stokes parameters of the incoming 
light. For matrix ( fTOb we find 



Qmrr = cos(2p)g + sin(2p)U and 
U corr = cos(2p)U - sin(2p)Q . 



(32) 



That is, a rotation of the angle reference as the one caused by the 
HWP offset cannot be corrected at all, while for a metallic reflec- 
tion U is still affected by a factor that depends on the material 
constants. 

Furthermore this method can only correct instrumental ef- 
fects for the part of the instrument that is rotated by 90°. Because 
in our case the rotated part is the analyzer (a rotation of the HWP 
is equivalent to a rotation of the analyzer), this procedure only 
corrects for the IP of the light train after the HWP. Another pos- 
sibility is to rotate NACO with the rotation adapter as a whole, 
as was done for a number of observations in the ESO archive. In 
both cases the IP of M3 remains uncorrected. Nevertheless, us- 
ing the "channel switch" method, the effects of the analyzer and 
the flat-fielding, which are difficult to quantify as described in 
section l2~5l and even parts of the IP of NAOS can be eliminated. 
We propose the following strategy. 

To archive the best attainable accuracy we propose to realize 
the 45°-switching between Q and U with the HWP, and the 90°- 
s witching with a rotation of the entire instrument (NACO). We 
now consider only the effects after M3, i.e. the system Tr x 
T(-/3) x 7(90°) xM NA0S x T(-90°). Then the numerical analysis 
shows that Eq.[29l gives 



Q corr = 0.894g- 0.387 U 
U cor r = 0.447Q + 0.775 U . 



The situation gets much easier if we measure at offset corrected 
angles (without HWP offset, the angles can be obtained from 
Table [TJ. Then we get 



Qcorr — Q 

U corr = 0.865 U . 



(34) 



Here we just have to additionally correct for the factor in U and 
the IP of M3, which reduces the number of free parameters of the 
model to 3. The effects of all the parts after M3 can be reduced 
like this to one factor. 

This method in particular eliminates the mentioned uncer- 
tainties of the parameters of the analyzer. Since the description 
of the parallactic angle-dependent part caused by M3 proved to 
be very accurate (a few tenth of a percent in polarization de- 
gree) in comparison to the remaining systematic deviations be- 
tween model and standards (about 1% in polarization degree), 
we expect that the proposed method will allow us to improve 
the accuracy by about a factor ten. We will investigate this in 
future calibration runs. We emphasize here that the "channel 
switch" method needs more than twice the time to obtain one 
set of Stokes parameters and is not suited for Sgr A* and its fast 
variability with time scales down to a few minutes. Another dis- 
advantage is the restriction of the field-of-view, which is caused 
by the rotation of NACO with the rotation adapter. 

4.2. The "boot strapping" method 

For time-resolved measurements in the crowded Galactic cen- 
ter field we have commonly used a "boo t strapping" metho d 



to calibrate the polariza tion data (e.g. lEckart etaL l2006allcl 



(33) 



Zam aninasab et al]|2010l) . This me thod has also been success- 
fully applied to wire-grid data (e.g. lEckart et alJ[T995l lOtt et al.l 
1999). In the presence of crowded fields with many weakly and 
only a few strongly polarized sources, it has the advantage of be- 
ing applicable without the availability of extensive data on cal- 
ibrator sources. Below we describe this method and investigate 
its uncertainties in detail. 

For the polarimetric "boot strapping" calibration of the light 
curves of Sgr A*, each channel is flux-density-calibrated with 
reference stars in a region of 2 arcseconds diameter surround- 
ing Sgr A* assuming total intensity brightness for each star. The 
sums of the orthogonal channels for 0° and 45° are averaged and 
taken as total intensity. With this total intensity and the galac- 
tic foreground polarization of 4%@25° dKnacke & Cappsll 19771) 
one obtains with Eq. [3] the expected flux densities for each star 
and channel. These flux densities are then compared with the 
time-averaged fluxes of the light curves of each star and channel, 
and a correction factor for each channel is obtained by averag- 
ing over all stars. Following this procedure the stars in the near 
surrounding of the GC show in average the foreground polar- 
ization and every source with similar polarization is calibrated. 
The value of 4%@25° is an average for the s ources toward the 
centra l arcsecond that has been measured by Knacke & Capps 
(11977b with arcsecond resolution. We could confirm this mea- 
surement by our independent calibration, which results in an av- 
erage of 5.3%@27° for the central arcsecond (Table©, which is 
equivalent within the errors. 

By using this calibration procedure one assumes that the IP 
affects the measurement by introducing different weighting fac- 
tors to the flux measurements of the four polarimetric channels, 
very similar to the assumption in Eq.|28] Here these factors are 
considered to be independent of the polarization of the consid- 
ered source (whereas the factors for the "channel switching" are 



11 



G. Witzel et al.: The instrumental polarization of NACO 



assumed to be independent of an instrument rotation). In partic- 
ular it cannot correct for an angular offset like the one caused 
by the HWP if the calibrator's polarization angle is significantly 
different from the polarization angle of the source one aims to 
calibrate (this offset corresponds to a Q <-> U cross talk, and 
the correction factors for each channel in this case depend on 
the direction of the linear polarization). We intend to answer the 
question of how this systematic error of the "boot strapping" cal- 
ibration influences the variability of the polarimetric light curves 
ofSgrA*. 

First we investigate the syste matic errors for a t heo retical 
light curve pattern as deduced in Zam aninasab et al] d2010l) for 
a polarized orb iting hot spot in an accret ion disk around Sgr A*. 
In Fig. 37 of IZamaninasab et all (120101) an apparent view of a 
hot spot in a Keplerian orbit at the innermost stable circular or- 
bit (ISCO) of a spinning black hole (with spin parameter of 0.5) 
is shown. The model predicts that the observer witnesses a mag- 
nification in flux according to lensing and boosting effects. The 
polarization angle on the observer's sky sweeps shortly before 
the total flux reaches its maximum, while the degree of polar- 
ization follows this maximum (Fig. fTOb . The existence this a 
pattern is an indicator for the strong gravitational regime and 
can be used as a tool for measuring the spin of a black hole 
dZamaninasab et all l2010h . While this model predicts the de- 
scribed pattern as a function of the normalized orbital time scale, 
we here set the orbital time scale to 30 min (as observations of 
Sgr A* suggest) and the center of the pattern (here at t=15 min) 
to hour angle zero, where the variability of the IP is strongest. We 
compute the Stokes vector for a source that shows the apparent 
polarization through the foreground of 5.3%@27° (see Tabled, 
apply the Mueller matrix for NACO to this vector, and deduce 
the normalized flux in each channel. These fluxes are then com- 
pared with the expected fluxes for a source o f 4% @ 25° (without 
IP) as assumed in previou s publications (e.g. Zamaninasa b et al.l 
2010. lMever et al.l l2006b). For each time the obtained correction 
factors are applied to the channel fluxes that have been calculated 
from the theoretical polarization pattern; this pattern had to be 
transformed before with M^aco to describe the actual measure- 
ment at the detector as predicted by our IP model. 

The resulting light curves are shown in Fig.[l0] The peak val- 
ues of the polarization degree are systematically underestimated 
by about 10% at an expected peak of 70%. Typical deviations 
are on the order of <5% for the degree of polarization and < 13° 
for the polarization angle. However, the angle mainly shows the 
expected HWP offset of about 13°, while other polarization ef- 
fects of the instrument have much smaller influences. Indeed, the 
resulting light curves look very much the same if the simulation 
only takes the HWP offset into account. For very low polariza- 
tion degrees the angle is ill defined. At these states real data with 
white noise contribution do not allow for detecting significant 
polarized flux, and therefore the interpretation of the polariza- 
tion angle is not possible in either way. The overall behavior of 
the variability is conserved. A compensation for the HWP offset 
during observations would eliminate almost all the effects intro- 
duced by this calibration method. 

4.3. Effects on time-resolved polarimetric measurements of 
Sgr A* 

As a final step we investigate the difference of the calibra- 
tion methods with respect to the observed polarimetric data of 
Sgr A*. For comparison a proper error estimate of the individual 
data points is needed. For the common "boot strapping" method 
the error was deduced from ( 1 ) the statistical variations of a com- 




5 10 15 20 25 30 



time [min] 



o 




1 5 10 15 20 25 



time [min] 

Fig. 10. Simulation of systematic calibration artefacts in light 
curves. The green curves show a typical p attern for polariza- 
tion d egree and angle as a function of time (IZamaninasab et ail 
120101) . The blue curves show the same data calibrated by the de- 
scribed "boot strapping" method. The difference between both is 
displayed in red. 



parison star near Sgr A* for each channel after flux density cal- 
ibration, and (2) the standard deviation of the correction factors 
for the calibration stars. This was performed by Gaussian error 
propagation for p and <p. Since the new calibration does not in- 
clude a flux density calibration, one has to estimate the statistical 
error of the photometry from the ADU (analog to digital con- 
version units) counts of a comparison star. Here it is important 
to first eliminate the correlated fluctuations of both orthogonal 
channels. This is achieved by subtracting one channel from the 
other after scaling the subtracted channel in a way that the aver- 
ages of both channels are the samfl The standard deviation of 
this difference is a good error estimate for the difference in flux 
between both channels (and the total flux as well), and can be 
propagated again. 

As a result we obtain Figures QT| IB. 21 and IB. 31 In general 
both calibration methods show very similar results within the 
statistical uncertainties of the measurements. The new calibra- 
tion shows a trend toward smaller polarization degrees, and the 
polarization angle shows a small systematic offset as discussed. 
Generally, the polarization angle is not well defined for small de- 
grees of polarization. The comparisons in Figures QT| IB. 21 and 
IB. 3 1 show that for total Ks-band intensities above 4 mJy and po- 
larized fluxes above 1 mJy the results of both polarization cal- 



7 For comparison stars of low polarization the scaling is not crucial. 
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Fig. 11. Total intensity (upper panel), polarization angle (mid- 
dle panel), and polarization degree (lower panel) of Sgr A* as 
a function of time (2006 June 1). The black data points are the 
result of the new calibration method, the red points are obtained 
by the previously used "b oot strapping" method as published in 
Zamaninasab et al.l (l2010h . The blue points show polarized flux 
(p x /). The meridian transit occurred at 32.4 min, 50.4 min be- 
fore the polarimetric measurements started. 



ibration methods are virtually identical. Only states of Sgr A* 
that agree with these conditions hav e been interpreted in the 
framework of a relativistic mod eling (Zamaninasab et al.ll2010l 
lEckart et alJl2006dlMever et al.ll2006bllah . 



5. Conclusions 

We presented a detailed analysis of the polarization calibration 
of the ESO VLT NAOS/CONICA system in the Ks-band. Using 
the Stokes/Mueller formalism for metallic reflections we intro- 
duced a polarization model of the camera/telescope system that 
excellently agrees with the measurements obtained on calibrator 
sources and sources in the Galactic center. We can qualitatively 
and quantitatively reproduce the instrumental polarization and 
show that a polarization angle offset of (13.2 + 0.3)° has to be 
taken into account for NACO data observed before autumn 2009. 
The model presented here enables us to correct for the instru- 
mental polarization of NACO with an accuracy better than 1% 
in Ks-band. Our approach allows us to extend the description of 
the IP to other wavelengths easily. It can also be applied to other 
telescopes and Nasmyth focus instruments. Additional calibra- 
tion observations may allow for an even better accuracy. The IP 
in Ks-band amounts to a maximum of 4%. We showed that the 
"boot strapping" method, which can efficiently be used in the 
crowded Galactic center field, yields the same results within the 
statistical uncertainties in bright flare phases when compared to 
the more exact and elaborate polarization model. 

In summary our investigation shows that for sources with 
statistical errors in Q and U significantly smaller than 0.5 % the 
polarization can be measured with an accuracy of better than 1 % 
in polarization degree. The accuracy of the polarization angle in 
these cases is < 5° for polarization degrees > 4%. For weaker 
sources like Sgr A* the accuracy of the polarization measure- 
ments is dominated by the statistical errors. We find that for to- 
tal intensity fluxes above 4 mJy and polarized fluxes larger than 
1 mJy the previously obtained results on the polarized emission 
from Sgr A* are unaffected by our new findings. Finally we pro- 



posed a calibration strategy that probably will allow for observa- 
tions with an accuracy of a few tenths of a percent in the future. 
For this strategy a compensation for the HWP offset while ob- 
serving is essential. On the other hand without offset even almost 
all of the effects introduced by the "boot strapping" calibration 
method would be eliminated, too. 
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Appendix A: Correction matrices for the optical 
components of NACO 

In this appendix we give an analytic solution for the correction 
matrix Cnaco m a wa y that makes it easy to program it. This 
solution was obtained by the inverse of the matrices in Eq. [13] 
which we also give here. The inverse of the rotation matrix in 
Stokes space (see Eq. ITOb can be obtained by changing the sign 
of the parameter p. The inverse of the Mueller matrix for metallic 
reflection ([8j) is given by 



r + r_ 

r_ r + 

Ac A s 

-A, A c , 



(A.l) 



with variables analog to lA.7l With the inverted Mueller matrices 
we are able to give an analytic correction matrix that transforms 
the measured Stokes vector into the Stokes vector on the sky 
(with the orientation east of north) 



S = C NAC0 x S ' = M'l x M- x on x S ' , 



(A.2) 



with M var the position depending and M co „ the constant part of 
the Mueller matrix for NACO 



AC 



with 



• bs p as p — bc p 



ecp + fs p c p + t fcp - eSpCp + u 
es„c„ - fcp + v 
- gs P 



ec l + f s P c p - w 
p 



dc„ 



-ds„ - gc p 





_ \alu 
u .i 'V 

Afl/M 

\alu 



(A3) 



and 



MZ 



with 



T+k - T-l (cp + spj T-k - T + lcp -T+hp 
-T+l + T_k(cp + Sp\ —T-l + T + kcp T + ksp 



(c p - sp) + m 



—T + msp 
-T+nsp 



T + mcp -T + n 
T+ncp T+m ) 



(A.5) 



k = (rff + (r!' 7 ) 
/ = 2rfr sU 
m = (Aff - (Af f 
and with 



n = 2Af Af 
T ± = r 2 ± r\ 



a lit 



A alu = 



1 r ± 

cos(6 alu ) 
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The parameters are 



p 




parallactic angle 


a 




rotator angle 


P 
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alu 

r \\ 




0.962622 


alu 
A. 
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^alu 




176.028° 


...v/7 
Ml 




0.98272 


r sil 
1 ± 




0.98872 


8 sa 




165° 


-goldj 

II 
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rl 




0.98484/0.981932 (S13/27) 


goldJI 




0.979588/0.97865 (S13/27) 


soldJI 




0.979642/0.980538 (S13/27) 



(A.6) 



(A.7) 



(A.8) 
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S P 
C P 
a 
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= sin(2p) 

= cos(2p) 

= sin(2/?) 

= cos(2/3) 

= -r ah ' cos(2ff) 
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-rf ' sin(2a) 
r a l " cos(2or) 
-Af" cos(2a) 
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( 

= Af cos(2a) / = -( 



rf+Af" )cos(2a) 



rf ' + Af") sin(2a) 



These matrices are not normalized, because it would make 
them more difficult to read. This just effects total intensity, the 
normalized Stokes parameters remain unaffected. To switch off 
matrix Tr l after a flat-field correction as described in 12.51 set 
r\ = r 2 . 



(A.4) 
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Appendix B: Tables and figures 

In Fig. IB. II we show the confirmation of our model through 



2007 May 15 



Stokes-fits for IRS16NW, IRS16CC, and S67. In Fig. |B3] we 
show two more Sgr A* light curves analyzed with the more ex- 
act and elaborate method presented here and the 'boot strapping' 
method used before. The data used for this paper are summarized 
in TablelB~T1 
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Fig. B.2. Polarimetric light curve from 2007 May 15. Plots ana- 
log to Fig.Qj] The meridian transit occurred at 86.5 min. 
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Fig. B.3. Polarimetric light curve from 2007 May 17. Plots ana- 
log to Fig.Qj] The meridian transit occurred at 126.4 min. 
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Fig.B.l. Stokes-fits (Q blue, U green) for IRS16NW (upper Fig. B.4. Orientations and angles for Eq.|2]and[3] with i// 
panel), IRS16CC (middle), and S67 (lower panel). See Fig. [9] 
and Table g] 
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Table B.l. Observations in the Ks-band. 



Source 


Date 


# Frames 


Time Interval (UT) 


Project ID 


PI 


Setup 


GC & Sgr A* 


2006 June 1 


156 


06:38:39 - 


10:44:27 


077.B-0552(A) 


Eckart 


Woll. & HWP (0°/45°), S13 


2007 May 15 


116 


05:29:55 - 


08:31:48 


079.B-0084(A) 


Eckart 


Woll. & HWP (0°/45 o ), S13 




2007 May 17 


192 


04:42:14- 


09:34:40 


079.B-0084(A) 


Eckart 


Woll. & HWP (0745°, S13) 


GC&IRS16 


2009 May 18 


286 


04:37:55 - 


10:19:54 


083.B-0031(A) 


Eckart 


Woll. & HWP (0745°), S13 


IRS21 


2005 July 30 


18 


06:44:50 - 


07:01:32 


075.B-0093(B) 


Eckart 


Woll. & HWP (0730760°), S13 


RCra88 


2004 June 11 


10 


10:15:49- 


10:25:49 


073.B-0084(A) 


Genzel 


Woll. & HWP (0745°), S13 


WD 1344 


2007 April 08 / 09 


10 






079.D-0441(B) 


Israel 


Woll. & rotator (0745°), S27 


WD2039-202 


2007 August 18 / September 21 


15 






079.D-0444(A) 


Israel 


Woll. & rotator (0745°), S27 


HD109055 


2009 March 10 / 14 


38 






082.D-0 137(B) 


Israel 


Woll. & rotator (0745°), S27 



Notes. The column "setup" shows the facilities used and at which angles the data were obtained. 
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